Entanglement in a three spin system controlled by electric and magnetic fields 



Jakub Luczak and Bogdan R. Bulka 
Institute of Molecular Physics, Polish Academy of Sciences, 
ul. M. Smoluchowskiego 17, 60-179 Poznan, Poland 
(Dated: April 16, 2012) 

We study the effect of electric field and magnetic flux on spin entanglement in an artificial trian- 
gular molecule built of coherently coupled quantum dots. In a subspace of doublet states an explicit 
relation of concurrence with spin correlation functions and chirality is presented. The electric field 
modifies super-exchange correlations, shifts many-electron levels (the Stark effect) as well as changes 
spin correlations. For some specific orientation of the electric field one can observe monogamy, for 
which one of the spins is separated from two others. Moreover, the Stark effect manifests itself in a 
different spin entanglement for small and strong electric fields. A role of magnetic flux is opposite, 
it leads to circulation of spin supercurrents and spin delocalization. 

PACS numbers: 71.10.-w, 03.67.Bg, 73.21. La 



I. INTRODUCTION 

The last decade has seen a great interest in applica- 
tion of concepts from quantum information theory, for 
example entanglement [1], to condensed matter theory 
[2]. Since entanglement represents unique quantum cor- 
relations, the concept has been applied to exploration of 
phenomena in strongly correlated many-fermion systems 
in order to gain insight into the nature of quantum phase 
transitions. In this paper we show how entanglement is 
related with a spin correlation function and how it can be 
controlled by an external electric and magnetic field. We 
choose a system of three coherently coupled semiconduct- 
ing quantum dots with three electrons, because it can be 
viewed as a realization of a three qubit system, which has 
recently been of great interest [3]-[7]. In such the system 
one can find two classes of truly three-partite entangled 
states, represented by the Greenberger-Horne-Zeilinger 
(GHZ) and the Werner (W) states [Il[8]. Recently there 
are attempts to measure and control these states [9 , as 
well as to apply them in logical gates [10 . In this paper 
we follow a scheme for universal quantum computations, 
proposed by Di Vincenzo et al. [5 for a spin system 
with exchange interactions in quantum dots, in which 
logical qubits are encoded in the doublet subspace with 
= +1/2 (see also [6]). 

Recent experiments demonstrate that in three quan- 
tum dots one can perform coherent spin manipulations 
pTNT3] . It is well known that spin manipulation can be 
controlled by electric field, for example in the systems 
with the spin-orbit interactions [Ml [15] , by applying in- 
homogeneous static magnetic field [161 Ull: by a light- 
induced magnetic field through the dynamical Stark ef- 
fect [18 or Raman transitions [19] (see also [20]). In our 
approach a role of the electric field is different, it modifies 
super-exchange coupling. In the system with a triangu- 
lar geometry the electric field breaks its symmetry and 
changes the quantum correlations between spins. Role of 
a magnetic flux is different, it induces spin supercurrents 
flowing around the triangular ring and is the main de- 
coherence source in that kind of systems [21 . One can 



expect that the magnetic flux acts destructively on the 
entanglement. 

The paper is organized as follows. In Sec. II we show 
that the concurrence, as a measure of entanglement, has 
an explicit relation with spin correlation functions and 
chirality in the triple spin system. Therefore one can 
have a simple interpretation of separability, monogamy 
and dark spin states. Sec. Ill describes our system within 
the Hubbard model and its canonical transformation to 
the Heisenberg Hamiltonian. We show that the electric 
field breaks the symmetry of the system and modifies 
exchange coupling, whereas the magnetic flux generates 
spin chirality. Detail studies of the concurrence, as well as 
the spin correlation functions and spin chirality, are pre- 
sented in Sec. IV. For a special orientation of the electric 
field we find a biseparable state (monogamy). Sec. V 
summarizes the paper. 



II. SPIN CORRELATION FUNCTIONS AND 
MEASURE OF ENTANGLEMENT 

We begin defining wave functions for three electrons 
in the three qubit system. These functions can be con- 
structed by adding a third electron to the singlet or 
triplet state (see [22 ). Two spin subspaces are possible 
to define: quadruplets and doublets. The quadruplets 
are the states with the quantum spin number S = 3/2, 
Sz = {±3/2, ±1/2} and the corresponding wave func- 
tions are constructed from a triplet state in the form: 

lQ-1/2) = ;^(4t4i4i + 4i4t4i + 4i4i4t)i^)' 
IQ1/2) = ;^(4i4t4t ^ 4t4i4t ^ 4t4t4i)|o)' 

IQ3/2) = 4t4t4tio)- 

Symbols c]^ are creation operators of an electron with 
the spin a in the qubit i acting on the vacuum |0). A 
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linear combination \GHZ) = -^{\Q_s/2) + IQ+3/2)) is 
known as the GHZ state. Two other functions IQ+1/2) 
and \Q-i/2) are cahed W-states. These states are weh 
known in the hterature [T| l2l l8l l24l |25] . 

In this paper the studies are focused on the doublet 
subspace with the total spin S = 1/2. These states 
were proposed for exchange-interaction universal quan- 
tum computations [5]-[7]. In many cases, as the one con- 
sidered in the next part of the paper, the doublet is the 
ground state. We assume that the system is kept coher- 
ent for a time sufficiently long in order to perform an 
entanglement measurement, and we fix the z-component 
of the total spin Sz = +1/2 in further considerations. 
The wave function can be expressed as 

l^^D,/,) = ai\Dy2)l + ^^2| A/2)2, (1) 

where: 

1^1/2)1 = ^(4i4t - 4t4i)4tio)' (2) 
1^1/2)2 = ;^[(4;4t + 4t4i)4t - 24t4t4j|o)' (3) 

In general should be expanded including a linear 

combination with the states with double site occupation: 
4>[^4t4^l^^' "^^^ state 11^1/2)1 is prepared by adding 
third electron to the singlet state, whereas 11^1/2)2 can 
be prepared from triplet or singlet states [22 . This con- 
struction allows us for gaining insight into monogamy 
and biseparability [26l |27] . 

As a measure of the entanglement we take the concur- 
rence. In order to calculate it [28 , we define a reduced 
density matrix of a pair of electrons in the quantum dots 
i and j: Qij = Tr/e^, where i, j and k denote different 
quantum dots, and ^ is a density matrix g = \^d){^d\ 
for the doublet subspace. Next we derive a matrix: 

Rij = Qiji^y ^ ^y)Qij{^y ^ ^y)^ (4) 

where Gy is a Pauli matrix and the asterisk denotes com- 
plex conjugation of Qij. The concurrence is calculated 
as: 

dj max{0, Ai - A2 - A3 - A4}, (5) 

where Ai, A2, A3, A4 are square roots of eigenvalues of 
Rij in descending order. Cij can take values between 
zero (for separate states) and one (for fully quantumly 
entanglement states). Using this definition one can cal- 
culate the concurrence for the doublet representation ([T]): 

2 

C12 = -|q^2||V3q^i +0^2!, (6) 
2 

Ci3 = -|q^2||V3<^i - Q^2|, (7) 
C23 = ^|3a?-ai|. (8) 



Let us now calculate spin correlation functions in the 
doublet subspace ([l]): 

(Si • S2) = i[-\/3(aia; + a\a2) - 2\a2\% (9) 

(Si • S3) = ^[\/3(aia; + ajaa) - 2\a2W (10) 

(S2-S3) = J(-3|ai|2 + |a2|2). (11) 

In general the coefficients ai and 0^2 can be complex, for 
example in the presence of the magnetic fiux. We show 
later that in this case a spin supercurrent occurs with a 
nonzero value of chirality 

(Si-(S2 xS3))=i^(aia;-ata2). (12) 

Comparing these results with the concurrence Cij ([6|- 
(|8| one can find, after some algebra, the following rela- 
tion: 

C,, = ^|(S,-S,+iSi-(S2xS3))|. (13) 

This fact allows us to propose the expectation value of 
spin correlation functions and chirality as an alternative 
measure of the entanglement. 

In a multi-qubit system it is interesting to define 
monogamy of entanglement. If two qubits are fully 
quantumly entanglement then they cannot be correlated 
with the third one [29\ In this case one says about 
monogamy, which expresses the nonshareability of en- 
tanglement. For a three qubit system monogamy can be 
measured by the concurrence Ci(^jk^ between a qubit i 
and other qubits j, /c, which is given by the one-tangle 
^'iijk) ~ 4det^^, where Qi = Ti jj.g. For biseparable 
states C^i^.j^-^ = 0, which means that the qubits j and 
k are fully quantumly entanglement whereas the qubit 
i is separated. This quantity is related with the linear 
entropy SL^Qi) = 2[1 — Tr(^?)] = 4det^i. Monogamy of 
entanglement satisfies the relation 

C2,fc)>C2+Cii (14) 

proved by Coffman, Kundu and Wooters [28 . For our 
case, in the doublet subspace, we have the equality 
^iijk) ~ ^ij ^ ^ik^ could cxpcct for purc states. 

III. MODEL OF CORRELATED SPINS ON A 
TRIANGULAR SYSTEM OF QUANTUM DOTS 

In this section we would like to present a specific three 
qubit example, namely an artificial triangular molecule 
built of coherently coupled semiconducting quantum 
dots. We show first that an external electric field and 
a magnetic fiux can modify spin states, and later, in the 
next section, how these external fields infiuence the en- 
tanglement between spins. 
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Our system of three quantum dots (Fig{T]) is described 
by the Hubbard model 



H = ^{€i ^ Eer cos[0 ^ {i - l)27r/3]}ni^ 



(15) 



Here, = + E'er cos[^ + ~ l)27r/3] corresponds to a 
shift of a local single electron level in the electric field 
E. The polarization energy E • P = e E • ViTiia- = 
Eer cos [0 -\- {i — l)27r/3]n^cr, where: e - the electron 
charge, denotes a vector of the i-qubit position and 
is an angle between ri and E, n^^r - an electron number 
operator. Later for simplicity we put = and denote 
Qe = Eer. The second term in ([l5| describes electron 
hopping between nearest quantum dots in the presence of 
the magnetic fiux ^ enclosed in the triangle. According 
to the Peierls scaling an electron gains during hopping 
a phase shift (j) = 27T^/{hc/e). The Coulomb onsite in- 
teraction of electrons on the quantum dots is included in 
the last term. 




FIG. 1: Model of a system of Three Quantum Dots placed in 
an external electric filed E and a magnetic fiux 

Using this model we can calculate all properties of the 
system numerically. In particular for three electrons we 
calculate the spin correlation functions as well as the con- 
currence and infiuence of the electric field as well as the 
magnetic fiux. To understand the results we use a canon- 
ical transformation of the Hubbard Hamiltonian ( 15 ) 



to an effective Heisenberg Hamiltonian. Taking the hop- 
ping integral t and the electric field qe sls small parame- 
ters with respect to the Coulomb interaction one can 
get the effective Heisenberg Hamiltonian jjlj |32] : 



i<j 



Si • (S2 X S3). (16) 



The first term describes the superexchange coupling be- 
tween spins, for which the exchange parameter Jij can 
be calculated to the third order in t/U ^31j 

-A-/A-,i - A-M-,1 - A-,iA-,i) , (17) 



where Aij = {U + e j — e^). In the limit Qe one can 
get an explicit form of Jij\ 
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One can see that the second term is proportional to g'^ 
and corresponds to the quadratic Stark effect. When 
the electric field rotates, this term leads to oscillations 
with the period tt. The linear Stark effect is described 
by the third term, which corresponds to the period of 
oscillations equal to 2it. The linear term in J23 is always 
positive, whereas for J12 and J13 the linear terms can be 
negative. At ^ = one can see that J23 increases linearly 
with gE^ whereas the couplings J12 = J13 and they first 
decrease, and next increase quadratically for a larger gE- 
At = 4tcos0 we get J12 = J13 = J23 and the system 
becomes uniform once again. 



The second term in the effective Hamiltonian (16) de- 
scribes chirality of electrons in the presence of the mag- 
netic fiux. The term is connected with the Aharonov- 
Bohm effect and with the persistent currents moving 
around the fiux enclosed by the three quantum dot ring. 
The coupling parameter calculated to the third order in 
t/U is given by [31 



imtmj 
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im jm 
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ji jm 

• A"^A"M(21) 

which in the limit gE ^ U simplifies to the form = 
— 12t^ sin ^/[/^. This parameter depends on the electric 
field in higher order terms of the expansion, but we ne- 
glect them in our studies. 



Using the Heisenberg model ( 16 ) one can derive many 
physical quantities analytically, in particular energy for 
the quadruplet and doublet states presented in the second 
chapter. For \GHZ) and \W) state one has the energy: 



Eghz = Ew = 0. 



(22) 



These states are independent of the electric field because 
the spins cannot be transferred between quantum dots 
(due to the Pauli exclusion principle). In the doublet 
subspace [ 1^1/2)1, 1^1/2)2]^ the effective Hamiltonian is 
expressed as 



4 {J: 



12 ■ 



31 



4J23 + J31) Y(-^31- Jl2+iJx) 

- J12 - iJv) 



31J 



(23) 
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The eigenenergies are 
1 



2(^12 + J23 + Jl3)±^, 



(24) 



where A = [37^/4 + Jf^ + JI3 + Jf^ - J12 J13 - J12 J23 - 
-^13-^23]^^^- For the considered case the doublet states 
are below Eqhz and Ew The corresponding eigen- 
functions may be written as 



|A/2)l 



1^/2)2(25) 



where z± = (J12 -2J23 + Ji3±2A)/[V3(Ji3 - J12 -iJx) 



IV. INFLUENCE OF EXTERNAL FIELDS ON 
ENTANGLEMENT 

We calculated the spin-spin correlation functions (S^ 
Sj), the concurrence Cij 



and Ci(^jk^ for the Hubbard 
model as well as for the effective Heisenberg model 
(16). The results for both the models were the same, 
within the relative accuracy better than 0.01, for the pa- 
rameters used below. Therefore our analysis is performed 
for the Heisenberg model (16), which is much simpler. 



A. Role of electric field 

Let us first study how entanglement can be modified by 
an electric field only (in the absence of the magnetic field, 
i.e. for ^ = and = 0). The electric field breaks the 
symmetry of the system, induces polarization (altought 
its value is small due to the strong onsite repulsion U). 
The electric field also changes the concurrence Cy, Ci(^jk) 
and the spin correlation function (S^ • Sj). Figj2] shows 
these quantities for the ground state as a function of the 
angle 9 of the electric field with respect to the axes of the 
system. The left panels are plotted for a small electric 
filed {qe = 1), when the linear Stark effect dominates. 
The period of oscillation is 27r for Cij, Q(j7c) i^i'^j)- 
One can see in Fig|2^ that for ^ = the concurrence 
C12 = Cis = and C23 = 1. The spins in the qubit 2 
and 3 are fully entangled which means monogamy. The 
spin in the qubit 1 is separated from two others (see also 
the red solid curve Fig. 2b showing that Ci(23) = at 
= 0). The spin correlation functions (Figj2]3) (81-82) = 
(81 • 83) = and (82 • 83) = -0.75 which means that 
the spins 2 and 3 are in the singlet state. The wave 
function has the form: I^Di/s) = 1^1/2)1 • According 
to the classification [27 we have the simply biseparable 
states for which the density matrix can be written as 
p = Pi p23' At ^ = the electric field is oriented to the 
quantum dot 1 and the system has the mirror symmetry 
with the exchange couplings J12 = J13 < J23 [see (18)- 
(20)]. A similar situation occurs for l9 = 27r/3 and 47r73, 
only if gE < 4|t|. 
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FIG. 2: (color on-line) Dependence of the concurrence dj, 
Ci(^jk) and the spin correlation functions (S^ • Sj) vs. the 
angle of the electric field calculated for the ground state. 
The (red) solid curves represent: C12, C'i(23) and (Si -82); the 
(blue) dashed curves: C23, ^'2(13) and (S2 -83); the (black) 
dotted curves: C13, C3(i2) and (Si • S3). Left panel is plotted 
for a small electric field qe — 1, right one is for a large electric 
field qe — The other parameters are: = 20, = and t 
is taken as unity. 



For a large electric field qe > 4|t| the situation is dif- 
ferent and more complex (see the right panels in Figj2|. 
Now, the period of oscillations is changed due to the 
quadratic Stark effect [the second term in Eqs.(18)-([2Q|)]. 
For 6 = the ground state is 1^^1/2)2 and we do not ob- 
serve monogamy. The spin in the qubit 1 is separable 
from two others for ^0 close to 27r/3 and 47r/3 (see Ci(23) 
in Figj2^). The condition for monogamy is the mirror 
symmetry of the system. Using Eqs.([l8|-(20 ) and J12 = 

arccos(— 2t/^£;) for 



^13 < J23 one can get the angle Oq 

separability of the qubit 1. Figj2]F shows that at ^0 the 
quantum dots 2 and 3 are fully quantumly entanglement 
with (S2 • S3) = —3/4 (see the dashed curve). 



B. Role of magnetic flux 

In the presence of the magnetic flux (^ 7^ 0) chirality 
of the spin system [described by the last term in (16)] 
becomes relevant [32 -38 . Recently Hsieh et at. [38] pro- 
posed to use chirality for quantum computations. Their 
quantum circuits are based on qubits encoded in chirality 
of electron spin complexes in systems of triangular quan- 
tum dot molecules. The magnetic flux removes degen- 
eracy between the states with different orbital momenta 
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and leads to spin supercurrents circulating around the 
triangle [33l [34] . We take the expectation value of the 
operator [37] 



a = (4/V^)Si • (S2 X S3) 



(26) 



as a measure of chirality. [33] In the basis of the doublets 
([1]) and using (25) one gets the expectation value 



2A 



(27) 



in the ground state. It means that chirality depends on 
the splitting A between the doublets which can be con- 
trolled by the electric field as well. 

In the absence of the electric field (when all exchange 
couplings are equal Jij = J) the expectation value (Cz) = 
±1, because the supercurrent is oc sin0 and circulates 
clockwise or anticlockwise (for O<^<7ror7r<0<27r, 
respectively). For this case the spins are delocalized and 
the expectation values of the spin correlation functions 
(Si-Sj) = —1/4. The concurrence Cij = 1/3 and Ca^jk) = 
V^/3, which describes maximal spin mixing. 

In a general situation, in the presence of both the fields, 
one can expect a competition between localization and 
delocalization of spins. These processes should manifest 
themselves in the spin correlation function (S^ • Sj) (in 
Cij, Ci(^jk^) as well as in the expectation value of chirality 
{Cz)' Fig. 3 and 4 present plots for (Si • S2), {Cz) and 
Ci(23) in the 0-(j) plane for a small and large electric field, 
respectively. Close to (/) = and tt the spin supercur- 
rent (oc sin (j)) is small and the Stark effect dominates. In 
the other regions the magnetic fiux becomes relevant and 
{Cz) reaches its extremal value ±1 (see the plateau with 
small oscillations caused by the electric field in Fig. 3b). 
According to ( [27| ) {Cz) is inversely proportional to the 
energy gap A between the doublet states and reaches 
maximal values at the symmetry points = 0, 27r/3 
and 47r/3 (see also Fig. 4b). The functions (Si • S2) and 
C'i(23) are very sensitive to symmetry breaking caused 
by the electric field. They show changes when the elec- 
tric field becomes larger qe > 4|t|cos0 (compare Fig. 3 
and 4). A detail analysis of the plots shows two contri- 
butions to the Stark effect: a linear and quadratic ones. 
The magnetic fiux reduces the linear component, and the 
quadratic Stark effect dominates [see also Eqs.(18)-(20) 
for the exchange couplings Jij]. If increases, the spins 
become more localized and the amplitude of the Stark 
oscillations, seen in (Si • S2) and Ci(23), increases. The 
localization process is monotonic, we could not observe 
any drastic changes - in contrast to the situation in the 
region close to ^ = and tt when the spin correlation 
functions and the concurrence change drastically their 
characteristics in large electric fields. 



V. CONCLUSIONS 

Summarizing, we showed that entanglement can be 
controlled by the electric field in three spin system of 




FIG. 3: Plots of the spin correlation function (Si • S2), the 
chirality {Cz) and the concurrence Ci(23) in the plane and 
= 27r$/(/ic/e) for a small electric field qe = 2 and the other 
parameters the same as in Fig. 2. 



three coherently coupled quantum dots with a triangular 
geometry. The studies were focussed on bipartite entan- 
glement in the subspace of doublets with Sz = +1/2, for 
which the concurrence was related to the spin correlation 
function and the spin chirality, Eq.(13). This relation 



was exemplified for the Hubbard model and its canonical 
transformation to the effective Heisenberg model. The 
super-exchange coupling exhibits a linear and quadratic 
dependence on the electric field (the linear and quadratic 
spin Stark effect), which manifests itself in different peri- 
ods of oscillations of the concurrence and the spin corre- 
lation functions when the electric field changes its orien- 
tation. The competition between these two Stark effects 
leads to different characteristics of the concurrence for a 
small and large electric field. For a special field orienta- 
tion we found a biseparable state, for which one of the 
spins is separated from two others fully quantumly en- 
tangled (monogamy). For small fields one should direct 
the field precisely toward the quantum dot {Oq = 0, 27r/3 
and 47r/3) to get the spin separation at this dot. In the 



6 



gE=6 




0/271 1 */27i 



FIG. 4: Plots of (Si • S2), (a) and Ci(23) in the O-cf) plane 
for a large electric field qe = 

case of the large electric field its orientation Oq should 
be different: the angle Oq depends on the strength of the 
field and it should be directed toward one of the opposite 
quantum dots. 

We also considered a role of spin chirality on entangle- 



ment. The magnetic flux ^ induces circulation of spin 
supercurrents and leads to spin delocalization. The bi- 
partite concurrence becomes uniform. For small electric 
fields the spins are delocalized, the Stark effect can be 
only seen in a very narrow range of ^. For larger fields 
the Stark effect becomes more visible, the concurrence 
and the spin correlations exhibit oscillations as a func- 
tion of the angle 9 of the electric field. Analyzing the 
oscillations one can see how the magnetic flux modifies 
relative contributions of the linear and quadratic Stark 
effect to entanglement. 

Our studies can be also related to recent experiments 
on coherent spin manipulation in three quantum dots 
pJT[T3] . We showed that the scheme proposed by Di 
Vincenzo et al. [5 , with logical qubits encoded in the 
doublet subspace 1 1^1/2)1 and 1/^1/2)2, can be realized 
due to the spin Stark effect, in which the electric field 
changes spin entanglement. The ground state is a super- 
position of the doublet states, which can be controlled 
changing orientation of the electric field. Therefore, this 
effect can be used to preparation of a proper initial quan- 
tum state and control a logical operation. Let us point 
out that Di Vincenzo's et al. scheme is different than 
operations between quadruplet and doublet states per- 
formed by Gaudreau et al. [13] in the experiment in triple 
quantum dots and between the singlet and triplet states 
in double quantum dots [39], where passages were ac- 
companied by reorientation of nuclear spins in quantum 
dots. Within Di Vincenzo's et al. scheme the dynamical 
passages are performed only between the doublet states 
(without nuclear spins) and total spin = 1/2 as well as 
its z-component Sz = 1/2 are conserved. 
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